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1. Introduction 

The B physics programme continues to play a crucial role in testing the CKM mechanism of 
quark flavour mixing and in determining fundamental Standard Model parameters. In recent years 
accurate measurements of numerous B physics observables have been realized, which calls for an 
equal improvement in refining the theoretical expectations within and beyond the Standard Model. 

The main obstacle for precise theoretical predictions are the complicated strong-interaction 
effects encoded in the hadronic matrix elements. The QCD dynamics often simplifies consider- 
ably in the heavy quark limit mj, S> Agco, which allows to establish factorization theorems that 
disentangle short- and long-distance effects. This separation provides the key for a systematic im- 
provement of the theoretical predictions by computing higher order radiative corrections, which 
should be supplemented by similar progress in the determination of the remnant non-perturbative 
hadronic parameters. 

Here we report on recent progress in the perturbative calculations for charmless hadronic and 
semileptonic B meson decays. As a technical account of these calculations has already been pre- 
sented in [jl|], we focus here on the phenomenological implications of the NNLO corrections. 

2. Hadronic B decays 

Most of the observables at current and future B physics experiments are related to hadronic 
two-body decays. Among these time-dependent and direct CP-asymmetries in penguin-dominated 
decay modes are of particular phenomenological interest due to their sensitivity to New Physics. 

The control of the strong-interaction dynamics in non-leptonic decays is obviously demanding. 
The factorization formula for the hadronic matrix elements of the operators in the weak effective 
Hamiltonian takes a twofold structure 



which consists of universal non-perturbative parameters (form factor F BM (q 2 = 0), decay constants 
/m, light-cone distribution amplitudes 0m) and perturbative hard-scattering kernels Tf' 1 , that con- 
tain the short-distance dynamics of the flavour-changing quark transition. The latter are currently 
being worked out to NNLO, i.e. at ff(aj). While the full set of hard-scattering kernels from spec- 
tator scattering (Tf 1 ) for tree ^ and penguin amplitudes [Q] is now available at NNLO, the ones 
related to the vertex corrections (Tf) are known to date for the tree amplitudes only \^ g]. 

The NNLO calculation is particularly important for direct CP asymmetries that are first gen- 
erated at &(cc s ). As in any perturbative calculation, it may thus help to reduce scale ambiguities 
of the leading contribution. The NNLO corrections may even change the pattern of CP asymme- 
tries significantly, since they can potentially be enhanced by large Wilson coefficients (which is not 
possible at even higher orders since the NNLO terms already have the full complexity). 

The current status of the NNLO calculation does not yet allow to discuss CP asymmetries. We 
may, however, already consider (CP-averaged) branching ratios of tree-dominated decay modes 
that do not depend significantly on the penguin amplitudes. As these observables are likely to 
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be dominated by their Standard Model contributions, they may serve as important probes for our 
understanding of the strong-interaction dynamics in non-leptonic decays. 

The NNLO analysis of the eleven tree-dominated B — » Tin /Tip / pp decay modes has been 
presented in [^] (for a similar analysis cf. In general colour-allowed decay modes turn 

out to be under much better theoretical control than colour-suppressed modes. Let us illustrate 
this point at the amplitude level: for the colour-allowed amplitude in the Tin channels one finds 
a\(KK) = 1.013^;^ + (+0.027+°;^)/, which is to be compared with a 2 {nn) = 0.195+^ + 
(— 0. lOl^Q Qgj)/ for the colour-suppressed one. It is striking that the latter suffers from substan- 
tial theoretical uncertainties, which can be traced back to a strong (and unfortunate) cancellation 
between different terms in the perturbative expansion. This makes the real part of a 2 particu- 
larly sensitive to the spectator scattering mechanism, which is normalized by the hadronic ratio 
fnfB/^BFf^iO). The poor knowledge of the B meson parameter I/Xb = J^dco/co 0b (ft)), in par- 
ticular, makes the theoretical prediction of the colour-suppressed amplitude rather uncertain. 

In order to test the QCD dynamics in hadronic decays it is useful to consider ratios of de- 
cay rates rather than absolute branching fractions. Particularly suited are ratios that involve the 
differential semileptonic decay rate at maximum recoil, 

^m 3 {MiM 2 ) = \ • (2.2) 

dF(B° -> Mji~ Vi)/dq 2 \ q 2 =0 

Experimentally this requires to measure the semileptonic decay spectrum over a sufficiently large 
number of g 2 -bins. Assuming specific parameterizations for the form factor shapes, the spectrum 
may then be extrapolated to q 2 = 0. At present this information is available for B — > n£v decays 
whereas the data on the B — > plv spectrum does not yet allow for an accurate extrapolation. 

In Table [T] we confront the NNLO prediction of the & n -ratios with experimental data. We 
stress that the theoretical predictions from Table [l] are based on a default set of hadronic input 
parameters (specified in Table I of |7|]), that is motivated by recent lattice and sum rule calculations. 
We see that the theoretical predictions are in good agreement with the data, which strongly supports 
the factorization assumption 1 . This is in particular true for the ratio (7T~7r ), which does not 
depend on the QCD penguin amplitude and on weak annihilation contributions at all; it thus gives 
clean access to \ai(7i7i) + Gt2(7i7i)\ 2 [Bp. Taking current data at face value, we may conclude that 
the colour-suppressed amplitude is somewhat enhanced, which may hint at a lower value of the B 
meson parameter X B c± 250MeV (the default choice adopted in [0] is X B = (400 ± 150)MeV). It 
would be interesting to verify if this conclusion is supported by the according ratio in the p -sector. 











^n{n+p-) 


Theory 
Experiment 


0.7Qj£g 
0.81^ 


o.8o^ : l| 


1 71+ - 27 
1 57+ 032 

1 - J/ -0.32 


9 77 +0.32 
■ L -' ' -0.31 

^■^-0.47 



Table 1: Ratios of hadronic and differential semileptonic decay rates in units of GeV . 



'The agreement is less pronounced for the ratio ,'% k {tz + tz^), which shows a much stronger dependence on the QCD 
penguin amplitude and hence on the specific input value for the weak phase y. This ratio is thus not particularly suited 
to test the dynamics of the tree amplitudes. 
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As long as the experimental information on the semileptonic B — > p£v spectrum is absent, we may 
instead consider ratios of two hadronic decay rates, 

R(MiM 2 M3M4) = — L ^r. (2.3) 

The ratio R(p L pl/p£p L ) ~ \a x {p L p L ) + cc 2 {PlPl)\ 2 /2\ai(p L p L )\ 2 yields complementary infor- 
mation on the tree amplitudes from the p -sector . One should keep in mind, however, that this ratio 
receives corrections from the QCD penguin amplitude and from weak annihilation in contrast to the 
semileptonic ratio £% p (p£~p£). From the numbers in Table ^| we infer that the NNLO prediction is 
again found to be smaller than the experimental value, which supports the hypothesis of enhanced 
colour-suppressed amplitudes and hence a lower value of Ag. 

Let us finally comment on the colour-suppressed modes, which are more complicated due to 
their strong dependence on hadronic input parameters. In contrast to the colour-allowed modes, it is 
in particular not possible to reduce these uncertainties by considering semileptonic ^ M -ratios since 
their dependence on |V M /,| 2 |F+ M (0)| 2 is weak. One may instead try to resolve the correlation among 
the theoretical uncertainties by considering hadronic ratios of two colour-suppressed modes. As can 
be seen in Table ||, this does unfortunately not lead to an improvement for the ratio that involves the 
7r 7T decay mode 3 . This is different for the ratio R(n°p° / p|'Pl)> which is less contaminated by the 
QCD penguin amplitudes. Consequently, the dependence on Xb and, somewhat accidentally, the 
one from the modelled power corrections drop out to a large extent. A more precise experimental 
value for this ratio may therefore give further insight into the role of power corrections in non- 
leptonic decays. The dynamics of the colour-suppressed amplitudes, however, should be probed 
with the cleaner ratios M % {tTtP), & p (p L pl), & p (7r-p°) and &„(nPp-). 



3. Semileptonic B decays 

Semileptonic b — » u decays provide a measure of the CKM matrix element \V u b\. The current 
discrepancy between inclusive and exclusive determinations calls for further progress on both sides. 
Here we report on the NNLO calculation for inclusive B — > X u £v decays in the so-called BLNP 
approach [|To|]. 

The theoretical description of inclusive B — > X u £v decays is complicated by the fact that exper- 
imental measurements have to introduce kinematical cuts to suppress the B — > X c £v background. 
This restricts the experimental information to the shape-function region in which the hadronic final 





R(PlP°l/P^Pl) 


R(K°p°/K°7i Q ) R{k {) p°/pIpI) 


Theory 
Experiment 


0.65^ : i? 
0.89^ 


1 5o +1 ™ 1 17+ - 45 
1 ? q+0.36 j qn+1.45 



Table 2: Ratios of two hadronic decay rates. 



2 The subscript L refers to the longitudinal polarization. 

3 Let us emphasize that the agreement between the central values is accidental for this ratio, since the theoretical 
prediction for the individual branching ratios are quite below the experimental data (the numbers can be found in ||7||). 
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state has large energy Ex ~ my but moderate invariant mass p\ ~ rn.bA.QCD- In this region of phase 
space a factorization formula for the structure functions has been put forward [ |TI| ] 

Wt ~ Hi JdcoJ(p 2 m )S(co), (3.1) 

which contains an universal non-perturbative quantity, the shape function S, and two perturbatively 
calculable objects, hard coefficient functions //, and a jet function J, that encode the short-distance 
effects. The NNLO calculation of the latter is now complete. While the two-loop corrections to 



the jet function have been worked out in [|12|], the hard coefficient functions have recently been 
computed to NNLO by various groups [13]. The latter calculation required to match the flavour- 
changing V — A current from QCD onto soft-collinear effective theory JI4]]. It has further been 
generalized to the tensor current, which finds applications in electroweak penguin decays JIB]]. 
The numerical impact of the NNLO corrections on the inclusive determination of \V u b\ has 



recently been analyzed in []16|]. Starting from the two-loop expressions for the hard and jet func- 
tions, the authors implemented the renormalization group improvement and a specific model for 
the shape function. From their analysis of partial decay rates, they conclude that the NNLO correc- 
tions can be significant. This statement, however, depends on the choice of the (arbitrary) matching 
scale ^ ~ (ra/jAgcz)) 1 / 2 - For fa = 1.5 GeV, which was the default choice in the earlier BLNP anal- 
ysis [|10|], the NNLO corrections are found to be important; they typically lower partial decay rates 
by about 15 — 20% while at the same time reducing the perturbative uncertainties. 

In their determination of \V u b\ the authors combine the NNLO prediction of the leading term in 
the heavy quark expansion with known power corrections up to ^(l/m|). From the experimental 
information on partial decay rates, which are based on different types of experimental cuts (lepton 
energy Ei, hadronic invariant mass Mx, hadronic variable P + = Ex — \Px\), they deduce sample 
values for \V u b\. Some of their results are collected in Table ||, which illustrate that the central 
values are shifted significantly at NNLO. Concerning the error estimate, the first one reflects the 
experimental uncertainty, while the improvement on the perturbative uncertainty can be seen in the 
second one. One further infers from the last error that the numerical value of the Z?-quark mass, 
which enters certain moment constraints of the shape-function model, has a large impact on the 
determination of |V„/,|. Given that the NNLO calculation has increased the discrepancy with the 
exclusive determination, further theoretical progress on the treatment of shape function effects is 



highly desirable (cf., for instance, Q17[ ] for an alternative implementation that is not based on a 
specific model). 



Exp. 


Method 


\v ub \ [io- 3 ] 

NLO 


\Vub\ [IO" 3 ] 
NNLO 


BABAR 

BELLE 

BABAR 


Ei > 2.0 GeV 
M x < 1.7 GeV 
P+ < 0.66 GeV 


■j Q7 + 29+0.37+0.26 
—0.23—0.25 
a cc , n 04+O.22+O.21 

3 30 + 23 +0 - 27+0 - 25 


4.30 ±0.241°;^ 
3.87 ±0.26^;™ 

3.55*0.241^1 



Table 3: Values of \V U b\ deduced from different experimental measurements of partial 
B —> X u £v decay rates (the numbers are taken from [|l6||). 
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